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Relativistic heavy ion collisions have reached energies that enable the creation of a novel state 
of matter termed the quark-gluon plasma. Many observables point to a picture of the medium 
as rapidly equilibrating and expanding as a nearly inviscid fluid. In this article, we explore the 
evolution of experimental flow observables as a function of collision energy and attempt to reconcile 
the observed similarities across a broad energy regime in terms of the initial conditions and viscous 
hydrodynamics. If the initial spatial anisotropics are very similar for all collision energies from 
39 GeV to 2.76 TeV, we find that viscous hydrodynamics might be consistent with the level of 
agreement for V2 of unidentified hadrons as a function of pT- However, we predict a strong collision 
energy dependence for the proton V2{pt). The results presented in this paper highlight the need 
for more systematic studies and a re-evaluation of previously stated sensitivities to the early time 
dynamics and properties of the medium. 

PACS numbers: 25.75.Dw 



I. INTRODUCTION 



II. ENERGY DEPENDENCE OF PERFECTION 



A physics case has been made that collisions at the Rel- 
ativistic Heavy Ion Collider (RHIC) produce a strongly- 
coupled system that evolves as a nearly perfect fluid 
[H; that is, the medium has a value of shear viscos- 
ity to entropy density 77/s near a conjectured minimum 
bound . This categorization of the quark-gluon plasma 
as a nearly perfect fluid has opened up interesting connec- 
tions to other strongly-coupled systems in nature 
Many features of the RHIC experimental data are well- 
described by theoretical calculations using the assump- 
tion of strong coupling. The most compelling of these 
results are fully relativistic viscous hydrodynamic calcu- 
lations that describe the momentum anisotropy patterns 
measured by experiment. The initial overlap geometry 
in these nuclear collisions has a significant eccentricity 
for non-zero impact parameters. Both the general ini- 
tial anisotropy and the event-by-event variations due to 
fluctuations in the nucleon positions can be described in 
terms of various moments 62, £3, £4, ... 01 . If the created 
medium were a non-interacting gas of particles, these spa- 
tial anisotropics would have no mechanism to translate 
into the momentum distributions of partons and even- 
tually final state hadrons. However, in the limit of very 
strong interactions between the constituents (i.e. very 
short mean free paths), one expects substantial momen- 
tum anisotropics that might be describable via model- 
ing based on viscous hydrodynamics. These momentum 
anisotropies are often described in terms of measured 
Fourier moments of the azimuthal distribution of par- 
ticles - V2,V3,V4, ... ■ Extensive measurements of the 
even moments of these observables by RHIC experiments 
indicate a strongly fiowing fluid medium. Experimental 
efforts are now underway to study the odd Fourier flow 
coefficients that result from fluctuations, motivated by 
the important observation of Refs. 



The recently completed first LHC heavy ion run at 
y/SNN =2.76 TeV has provided an enormous increase in 
collision energy over the top RHIC center-of-mass energy 
of ^/SNN = 200 GeV. One of the more intriguing results 
is the measurement from the ALICE collaboration of the 
elliptic flow V2 as a function of transverse momentum 
for inclusive unidentified charged hadrons (ft,+ + h~) 
as shown in Figures [1] and [2] 
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FIG. 1: Experimentally measured «2{4} as a function of 
transverse momentum pr for unidentified charged hadrons 
in the 20-30% centrality selection. Published results from 
the ALICE experiment at 2.76 TeV ^ and the STAR ex- 
periment at 200 GeV [l3| are shown, along with preliminary 
results from lower energies [l2|. The lower panel shows all 
data sets divided by a common fourth-order polynomial fit to 
the ALICE data points. 



In comparing their measured W2{4} (a technic^uc for 
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FIG. 2: Experimentally measured W2{4} as a function of transverse momentum pr for unidentified charged hadrons for centrality 
selections 10-20%, 20-30%, 30-40%, and 40-50% from top to bottom. Published resuhs from the ALICE experiment at 2.76 
TeV [il and the STAR experiment at 200 GeV [H are shown. The right panels show the two data sets divided by a common 
fourth-order polynomial fit to the ALICE data points in each separate centrality selection. 
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measuring V2 via four-particle cumulants [l^), with pre- 
vious measurements at lower energies at RHic 
a striking agreement is seen. The optimal approach to 
quantify the level of agreement between the values of 
V2{pt) measured at LHC and RHIC would employ a full 
statistical and systematic uncertainty constraint fit fol- 
lowing the methods developed in Ref. [l^. However, we 
note that in fitting a fourth-order polynomial to the pub- 
lished ALICE data, which include only statistical uncer- 
tainties, we obtain a /d.o.f. = 3.4/13 corresponding 
to a p-value = 0.996; any inclusion of systematic errors 
would only decrease the and increase the p-value. The 
very low values of result from the fourth-order cumu- 
lant method which produces statistical correlations be- 
tween the data points [l^ . Additionally, other sources of 
systematic uncertainties are not fully quantified for either 
ALICE or STAR data sets. Regardless, as shown in the 
lower panel of Figure [T] and the right panels of Figure [5J 
the agreement between the 2.76 TeV and 200 GeV data is 
at the level of a few percent from px = 0.5-2.5 GeV, with 
larger deviations possible at lower and higher pT- For the 
20-30% centrality bin shown in Figure [1] the agreement 
persists at the 5% level down to energies at least as low 
at 39 GeV. 

The PHENIX experiment has published data in Au- 
Au collisions a t J snn ~ 200 GeV on V2 using the event 
plane method [l5| that reveals a similar agreement [l6| . 
However, caution is warranted regarding too fine a level 
of agreement or disagreement on account of the different 
sensitivities of the methods to fiow fluctuations, which 
are appreciable at the 10-20% level in the final measured 
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III. PERFECT OR TOO PERFECT 

In this paper, we discuss some of the implications and 
new questions raised by this agreement of the V2 measure- 
ments. In doing so, we consider the following scenario. 
Imagine that we could prepare a system with a specific 
spatial distribution of energy density (in particular with 
a specific set of spatial anisotropics) and immediately 
evolve the system as a fluid with a speciflc flxed value for 
rj/ s to some appropriately large time. Let the medium 
have an Equation of State (EOS) where the speed of 
sound is independent of temperature - for example an 
ideal non-interacting gas with e = p/3. After the long 
evolution time, the system immediately breaks up into 
hadrons via the Cooper- Frye freeze out mechanism p^ . 
One might expect that the momentum anisotropy (e.g. 
V2) as a function of pt for hadrons would be the same re- 
gardless of the initial energy density scale. To be speciflc, 
if the central energy density were a factor of four higher, 
but the spatial distribution of that larger energy density 
were identical, then despite larger pressure gradients in 
all outward directions, the relative pressure gradients in 
different directions would be the same. Thus, if one spec- 
ulates that the fluid created at RHIC is ideal {rj/s ~ 0) 



or nearly ideal {rj/s = I/Att) and the same is true at the 
LHC, perhaps this explains the near identical nature of 
the experimental data. A number of assumptions (some 
more justifled than others) went into the above scenario. 
We examine some of these assumptions within the con- 
text of various models and speculate on what conclusions 
might be drawn. Useful reviews of recent work on viscous 
hydrodynamics can be found in Refs. [l"9l - [2^ . 




FIG. 3: Viscous hydrodynamic results using an ideal gas EOS 
and rj/s — 0.001. Shown are the V2 for unidentified hadrons, 
pions, and protons as a function of pr- The solid (dashed) 
lines are for the Ti = 420(340) MeV, and the lower panel 
shows the ratio of V2 values from the two cases. 

First, we use the publicly available viscous hydrody- 
namic model of Romatschke and Luzum [2^, to test 
this simple scenario. We fix the initial conditions for 
Au-Au collisions with impact parameter 6 = 7 fm using 
Monte Carlo Glauber results for binary collision posi- 
tions [ill. We use an ideal gas EOS, ij/s = 0.001 (not 
exactly zero for numerical stability reasons and known 
to reproduce the results of ideal hydrodynamical cal- 
culations), and an isothermal freeze-out temperature of 
140 MeV followed by Cooper-Frye hadronization and res- 
onance decays. We model two cases, one with initial tem- 
perature T, = 340 MeV and a second with T; = 420 MeV, 
as estimates for initial RHIC and LHC temperatures re- 
spectively We use identical initial spatial distribu- 
tions for both cases (with just a rescaling of the energy 
density for the appropriate initial temperature), in order 
to allow us to isolate the other dynamical effects. 

The hydrodynamic results are shown in Figure [3] As 
expected, the pions (and thus the unidentified hadrons 
which are dominantly pions) have the same V2 within a 
few percent for pT > I GcV. However, the proton V2{pt) 
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FIG. 4: Viscous hydrodynamic results using an EOS determined from lattice QCD and rj/s = 0.001 (left panel) and 77/s = l/An 
(right panel). Shown are the V2 for unidentified hadrons, pious, and protons as a function of pr- The solid (dashed) lines are 
for the Ti — 420(340) MeV, and the lower panels show the ratio of V2 values from the two cases. 



pattern appears shifted out in px for the higher initial 
temperature case. This is most likely due to the larger 
radial boost mapping onto the momentum distribution 
of protons. One speculation is that the deviation at low 
Pt for pions may be from resonance decay contributions. 
However, we have checked the primary (without decay 
contribution) V2 for pions and similar differences between 
the two initial temperature cases is observed. 

Second, we extend this question to the case of the QCD 
EOS (using the lattice QCD inspired EOS from Ref. U^) 
and for rj/s = 0.001. The left panel in Figure |4] shows 
the results from this calculation again with two initial 
temperatures of 420 and 340 MeV. Despite the variation 
in the speed of sound Cs as a function of temperature due 
to the EOS, we find a similar qualitative result. For the 
range pr = 0.5 — 3.5 GeV the pions, and therefore the 
unidentified charged hadrons, have a V2 that is the same 
within approximately 5%. It is notable that this V2 signal 
is actually somewhat stronger for the lower T,; = 340 
MeV. However, the proton V2{pt) pattern again appears 
shifted out in pT for the higher initial temperature case. 

The right panel of Figure |3] shows the results for the 
same QCD EOS but now with ij/s = I/Att. In this case 
there is a noticeable difference at the highest px for all 
particles with the lower temperature case having a lower 
V2- This observation suggests that even a minimal vis- 
cosity plays a significant role in limiting the growth of 
V2 with Pt for > 3 GcV at RHIC energies, while the 
significantly higher energy densities at the LHC extend 
the regime where inertial forces dominate dissipation. We 
have run the Ti = 420 MeV case to higher hadron pT and 
find that the V2 values do start to saturate for pT > 6.5 
GeV at a level near V2 ~ 0.44. See Section ITVl for a dis- 
cussion of other effects that must also be considered in 



this momentum regime. It is possible a hint of this is 
seen in the experimental 200 GeV data in Figure [2] in all 
centrality selections; but current uncertainties preclude 
any strong conclusion. Again, the proton V2 pattern is 
shifted to higher pT, leading to a rising ratio for the lower 
to the higher temperature case at pr < 2.0 GeV. Thus, 
the measured V2 for protons at the LHC is much antici- 
pated. 

Note that we have used smooth Glauber initial condi- 
tions for fixed impact parameter b = 7 f m and no hadron 
re-scattering after freeze-out of the fluid in these calcula- 
tions. Therefore, these quantitative calculations should 
not be compared directly with the experimental data, but 
rather reveal the qualitative trends with changing initial 
temperature. 



IV. SENSITIVITY TO 77/s 

A crucial experimental question is the sensitivity of 
V2 {pt) to the value of 77/s. A corresponding critical ques- 
tion for theory is the assumption employed by nearly all 
hydrodynamic calculations to date that rj/s is indepen- 
dent of temperature. In this section we examine these 
issues. 

We have calculated the 1)2 (pr) for unidentified charged 
hadrons for RHIC temperatures with different variations 
in 77/s as shown in Figure [5l The black solid curve cor- 
responds to rj/s = f/47r, the blue solid curve to 77/s ~ 0, 
and the blue dashed curve to 77/s = 2/47r. Note that we 
have not re-normalized the initial entropy density of the 
medium to maintain a fixed final particle multiplicity or 
mean transverse momentum as done in other works - for 
example see Ref. [l^. For 10% variations in 7;/s (green 
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FIG. 5: Viscous hydrodynamic results using QCD EOS for 
Au-Au b=7 fm events, initial temperature T = 340 MeV, and 
with varying rj/s. Shown are the V2 for unidentified charged 
hadrons as a function of pr- The black curve is for rj/s = 
l/47r. The other colors correspond to changes of 10, 20, 40, 
100% in 77/s. The solid (dashed) lines correspond to a decrease 
(increase) in the ratio. The lower panel shows the rato of all 
curves to the l/47r case. 



FIG. 6: Viscous hydrodynamic results for 112 versus px for 
unidentified charged hadrons using a QCD EOS and initial 
temperature T = 420 MeV. We compare results with fixed 
rj/s = l/47r, 77/s = 2/47r, and two cases with a temperature 
dependent value where rj/s = 1/Att for T < 340 MeV and 
rj/s = 2/47r for T > 340 MeV or 77/s = 2/47r for T > 180 
MeV. The lower panel shows the rato of all curves to the 
constant 7?/s = l/47r case. 



curves), there is very little change in the V2 values (less 
than 5%). In fact, for these small 10-20% increases in 
77/s, V2 shows a slight increase for pT below ~ 2 GeV. 

It is clear that even a ±40% change in 77/s (red curves) 
results in only a 10% variation in V2- Furthermore, it is 
important to note that this 10% variation is nearly in- 
dependent of pt from 0-5 GeV. Many comparisons of 
experimental data and viscous hydrodynamics focus on 
the differences at higher px because the theoretical cal- 
culations with different values of 77/s visually split (as 
shown on a linear scale). However, this region in px is 
also sensitive to other effects such as the implementation 
of the dep arture 5f{p) from the equilibrium distribution 
foip) and path length dependent jet energy loss (see 
for example 129|). 

Given the relative insensitivity of V2{pt) to the pre- 
cise value of 77/s, it is important to identify analysis pro- 
cedures that maximize the sensitivity to this important 
quantity. Plotting the ratio of V2{pt) between experi- 
ment measurements and theoretical calculations reveals 
important discriminating power at low pT- It is at this 
low Pt that the experimental statistical uncertainties are 
smallest, and the systematic uncertainties (which require 
further attention to fully quantify) are percentage uncer- 
tainties and thus largely independent of pt-- 

In most hydrodynamic calculations the specific 77/s 



value which is used is often referred to as an average. 
However, it is not simple to define what type of average 
this quantity is, as it is weighted in some complicated 
fashion over the space and time evolution. More realisti- 
cally, it should be regarded as simply a fixed value used 
to obtain a first-order estimate. In this simple picture, 
the RHIC and LHC mediums would require values for 
77/s that differ significantly less than 40%. 

We can examine the sensitivity to possible 77/s varia- 
tions with temperature by a straightforward procedure. 
As a baseline, we have run with the QCD EOS and an 
initial temperature of 420 MeV and constant values for 
77/s = l/47r and also rj/s = 2/47r. For comparison we 
have implemented a simple step-function temperature de- 
pendence where for T < 340 MeV 77/s = l/47r and for 
T > 340 MeV 77/s = 2 x l/47r (i.e. 100% larger). We 
also consider a second case where the step-function oc- 
curs near the transition temperature Tc = 180 MeV, that 
is, 77/s = l/47r for T < and 7;/s = 2 x l/47r for T > Tc. 
The results are shown in Figure [6] 

For the first case, the results show essentially no dif- 
ference in V2{pt) over the px range shown. The T=340- 
420 MeV range is that explored in the early times for 
the LHC created medium, and then after cooling it ex- 
plores the same temperature range as that for the RHIC 
created medium. Notably, there is no measurable differ- 
ence in V2{pt) for a factor of two change in the 77/ s ratio 



6 



for the early time higher temperature medium. Even if 
the step-function in the value of ry/s occurs at the much 
lower temperature of Tc = 180 MeV, then only a 10% 
decrease in V2 at all px is seen. The modest values of 
these changes is striking, and pose serious challenges to 
more precise extraction of transport coefficients via such 
measurements. 

Recently a study exploring a family of four 
temperature-dependent 77/s cases has been per- 
formed [30j . with a change in 77/s at the transition 
temperature Tc = 180 MeV that includes a possible 
sharp rise in 77/s just above the transition (for example 
in the case labeled HQ the 77/s value rises to 10 x l/47r 
by T > 400 MeV). While these very large values for 77/s 
call into question the validity of the calculation in that 
parameter space [3l|, even in this extreme comparison, 
only a modest (less than 15%) change in the predicted 
V2 is found at LHC energies (specifically the so-labeled 
LH-LQ to LW-HQ cases). Qualitatively this is consistent 
with our finding where our very small (by comparison) 
factor of two change located at a higher temperature 
(T > 340 MeV) has very little impact. Thus, it appears 
premature to decide whether the higher temperature 
region explored at the LHC (before cooling and evolving 
over the same temperature range as RHIC) has the 
same or different 77/s (as predicted in some models - 
for example [s^l)- It is clear that more detailed studies 
of the temperature dependence are necessary; work in 
that direction is underway [23, H^l- A parallel critical 
area of research is recent work towards calculating these 
dynamical properties (for example 77/s as a function of 
temperature) for QCD on the lattice [33| . 



V. INITIAL CONDITIONS 

In all of the above discussion we have deliberately as- 
sumed that the initial spatial eccentricity and thermal- 
ization times do not vary between RHIC and LHC ener- 
gies in order to separately study the effects due to 77/s 
variations, the EOS and particle mass. Since ideal hy- 
drodynamics predicts that V2 should be proportional to 
the initial eccentricity, it is also important to investigate 
whether the initial spatial distribution at the LHC is 
the same as at RHIC. Another expected difference be- 
tween these two energy regimes is the equilibration time 
T. It is often stated that the hydrodynamic matching to 
experimental data at RHIC indicate rapid equilibration 
T < 1.0 — 2.0 fm/c. If in fact the final V2 pattern is 
sensitive to the equilibration time, and r is significantly 
smaller at the LHC than at RHIC, the near-identity of 
the V2{pt) data from RHIC and the LHC is deeply myste- 
rious. Of course, one solution to this mystery is that the 
input assumptions are wrong, as discussed in Ref. [23j . 
where it is argued that the purported sensitivity to rapid 
equilibration is incorrect 

In order to understand the potential range of variation 
of initial eccentricity, we have examined two of the cur- 



rently used phenomenological models for calculating the 
initial spatial eccentricity. First, we use a Monte Carlo 
Glauber framework 3J] to calculate the different eccen- 
tricities as a function of collision centrality percentile 
as determined by the number of participating nucleons 
(Npart)- The most appropriate basis of comparison to 
compare data sets are the centrality percentiles (e.g. 20- 
30% of the total inelastic cross section) being used by ex- 
periments; the eccentricity ratios show considerably more 
variation if comparisons are made at fixed participant 
number or impact parameter. 

As already mentioned in Section HI there are multi- 
ple techniques for measuring V2 that give results vary- 
ing by of order 20%. These variations are due to the 
different influences of flow fluctuations and non-flow ef- 
fects [l3l ■ If we assume that initial eccentricities e„ trans- 
late into momentum anisotropics Vn in individual events, 
then the most applicable for the 7^2 {4} measurements is 

r 2 1 

£2(4} = 2(e^) — (e*) , which is shown in the left 

panel of Figure [T] One observes an approximately 5% 
(10%) larger eccentricity for Au-Au at 200 (39) GeV col- 
lisions compared with Pb-Pb collisions at 2.76 TeV. This 
5% difference between 200 GeV and 2.76 TeV values was 
previously noted in 

There are three important contributors to these dif- 
ferences. First, the simple difference of atomic mass 
between Au (197) and Pb (208). Second, the nucleon- 
nuclcon inelastic cross sections are significantly increas- 
ing across the energy interval spanned by RHIC and the 
LHC, ranging from 34 ± 3 mb at RHIC 39 GeV, 42 ± 2 
mb at RHIC 200 GeV, and 64 ± 5 mb at LHC 2.76 TeV 
(although this last value is not yet experimentally final- 
ized). A larger cross section reduces the fluctuations in 
the number of participants, particularly in the periphery 
of the interaction zone, which in turn reduces the magni- 
tude of e2{4}, as seen in the left panel of Figure[71 Third, 
and related to the second, is that the larger total A-A in- 
elastic cross section at the LHC changes the mapping of 
centrality percentile to impact parameter range. 

In the right panel we show the £2 values with respect to 
the reaction plane. This shows a slight larger 8% (15%) 
differences between 200 (39) GeV and 2.76 TeV. When 
comparing the corresponding values of €2 with respect 
to the reaction plane and e2{4}, there is an increasing 
level of disagreement for more peripheral reactions, which 
is understood to result from non-Gaussian fluctuations 
that reflect the un derly ing Poisson fluctuations from dis- 
crete nucleons [3J-[3g. More central collisions quickly 
approach the limit of Gaussian fluctuations where the 
two estimates agree. However, for some hydrodynamic 
model calculations the role of (presumably real) fluctu- 
ations has been greatly reduced in spite of using Monte 
Carlo Glauber initial conditions by keeping the events 
fixed with respect to the reaction plane and averaging 
over many events to create smooth initial conditions as 
input. Note that this is exactly what is done for the hy- 
drodynamic model comparisons in [37j . where even for 
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FIG. 7: Monte Carlo Glauber results for Au-Au collisions at 39 and 200 GeV and Pb-Pb collisions at 2.76 TeV. The left panel 
shows £2 {4} with respect to the reaction plane, and the right panel shows the RMS £2 with respect to the participant plane, all 
as a function of collision centrality percentile (determined from the number of participating nucleons). The lower panels show 
the ratio of the RHIC to LHC eccentricity values. 



40-50% central events this results is a 10% underesti- 
mate in the average eccentricity. The effect is even larger 
(more than a 25% underestimate) for the 70-80% cen- 
trality. Caution is therefore warranted in making such 
comparisons with hydrodynamic models using this ini- 
tialization and data. A useful discussion of this point is 
made in Ref . (ssj , and may point to the need for running 
multiple hydrodynamic events on individually fluctuating 
initial conditions (39| . 

We have also made similar comparisons for the Color 
Glass Condensate (CGC) model motivated initial condi- 
tion calculation. This class of Monte Carlo calculations 
start with the same Glauber model described above to 
determine the participants in a collision. However, a 
spatial region with say 20 participating nucleons is no 
longer assumed to produce ten times the energy density 
as a region with 2 participating nucleons. This is due 
to the assumed saturation in the number of low-x glu- 
ons, whose wavefunctions overlap across the longitudinal 
extent of the nucleus. These calculations produce eccen- 
tricities that are larger than those from the pure Monte 
Carlo Glauber based on participating nucleons alone jiO] . 

Here we utilize the rc-BK option [4l[ for the CGC ini- 
tial conditions. Figure [8] shows €2 {4} with respect to the 
participant plane as a function of centrality percentile for 
Au-Au at 39 and 200 GeV and Pb-Pb at 2.76 TeV. One 
observes in the lower panel that the eccentricities agree 
within less than 2% for much of the centrality range (10- 
70%). At first this result seems counterintuitive since for 
higher collision energies, one should be probing lower- a; 



gluons, saturation effects should increase, and the ec- 
centricity should increase relative to the Monte Carlo 
Glauber case. One does see a slight hint of this in that 
the Glauber difference between 200 GeV and 2.76 TeV 
shown in the left panel of Figure [7] is 5% and in this case 
it is reduced to 2%. 

In order to test this picture further, we have calculated 
the mean eccentricity with respect to the participant 
plane for Pb-Pb collisions for a large range in collision 
energies (17 GeV - 200 TeV). Shown in the right panel 
of Figure [5] are the resulting eccentricities for the rc-BK 
CGC case (solid lines) and the pure Monte Carlo Glauber 
(dashed lines) as a function of the number of participat- 
ing nucleons. In the lower right panel we calculate for 
each collision energy the ratio of CGC to Glauber eccen- 
tricity. As expected, the larger the collision energy, the 
larger the modification in the number of gluons and thus 
the larger the increase in eccentricity for the CGC case. 
However, the relative increase in eccentricity is only « 
10% from 200 GeV to 2.76 TeV. Curiously, this increase 
in eccentricity from saturation effects is largely canceled 
by the decrease in the underlying Monte Carlo Glauber 
eccentricity and in the mapping of centrality percentiles 
for Au-Au and Pb-Pb, as shown in the left panel. 

One striking feature is that between 17 and 200 GeV 
the CGC predicted ~ 20% increase of eccentricity relative 
to the Monte Carlo Glauber changes only very slightly. 
In fact, the authors of the model note that for energies 
below 130 GeV the formalism breaks down since one is 
probing gluons with x > 0.01. This raises the interest- 
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FIG. 8: The left panel shows the rc-BK CGC calculation for £2 {4} with respect to the participant plane for Au-Au at 39 and 
200 GeV and Pb-Pb at 2.76 TeV as a function of collision centrality percentile. The lower left panel shows the ratio of LHC to 
RHIC eccentricities. The right panel shows the CGC and Monte Carlo Glauber mean eccentricities for Pb-Pb collisions over 
a broad range of collision energies. The ration of CGC eccentricities to those for the Monte Carlo Glauber are shown in the 
lower panel. 



ing possibility that below some energy, there should be a 
transition in initial conditions from the larger CGC cal- 
culated eccentricities to the lower Monte Carlo Glauber 
values. Data from the CERN-SPS heavy ion fixed target 
program have not led to a definitive answer regarding 
V2 at these energies due to different methods, central- 
ities, and varying baryon contributions to unidentified 
hadrons [i^. With additional data from the RHIC en- 
ergy scan, it will be very interesting to see how and if 
the V2 decreases as the center-of-mass is lowered (as pre- 
dicted by many due to viscous effects in the hadronic 
re-scattcring stage). The possible difference in initial ec- 
centricity will require careful study (since CGC effects 
are not expected to play a role at these energies), and 
may be masked by other effects with a larger dynamics 
range, such as the (also poorly known) modifications to 
the EOS at finite baryon chemical potential. 



VI. AVERAGES AND KNUDSEN NUMBER 
METHOD 

The above discussion regarding the sensitivity to vari- 
ous parameters and their temperature dependence raises 
the question of whether simple models with a single rj/s 
value, constant speed of sound, and a single characteris- 
tic temperature can go beyond simple dimensionful scal- 
ing arguments to both capture the key physics and pro- 
vide a basis for quantitative extraction of key dynamical 



parameters such as rj/s. For example, it has been pro- 
posed that one can extract 77/s from a fit to V2 (integrated 
over all pt) versus particle density using a simple scaling 
ansatz in the Knudsen number K = X/R [i^ (here A is 
the mean free path for momentum transport and ^ is a 
characteristic system size discussed below). Previously, 
we had investigated some of the underlying assumptions 
for this model and found important ambiguities in the 
parameterizations and sensitivities |44|. Despite these 
findings, the quantitative values from [43| are frequently 
quoted as reliable estimates for 77/s. 

It is of course possible to test this Knudsen scaling 
formalism by comparing it to the results from viscous 
hyd rody namic models where one knows the exact input 
rj/s [45|. Those authors show that a Knudsen-based scal- 
ing of the results with 77/s exists, but did not check the 
formalism for self-consistency by using the same scal- 
ing [4^ to compare the implied value for rj/s to the input 
value used in the hydrodynamic calculations. We have 
done this following the exact procedure detailed in [i^l 
whereby one arrives at a compact relation: 



- = Q.32K0TR 

s 



Me) 



ih 



iv2/e) 



- 1 



(1) 



where Ko is a constant, T is the single constant tem- 
perature, R = {x^) + 1/ (y^) is the characteristic 
scale for the strongest gradient in the initial matter con- 
figuration, and {v2/e)ih is the ratio achieved in the ideal 
hydrodynamic limit. 
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FIG. 9: The left panel shows results for V2/^part as a function of the transverse particle density {1/ SrdN / dy) from calculations 
of viscous hydrodynamics with fixed rj/s — 0.001, l/47r, and 2/Att. The red and blue curves are fits to the calculation using the 
Knudsen number formalism. The right panel shows the resulting Knudsen number extraction for 77/s as a function of the same 
transverse particle density. The dashed lines represent the true rjjs for the two cases as input to the viscous hydrodynamic 
model. 



Shown in Figure [S] (left panel) are the results from 
viscous hydrodynamic calculations (using the model 
from [i^l) cast in these quantities for three differ- 
ent fixed (centrality independent) values of 77/s = 
0.001, l/47r, 2/47r. The two finite viscosity cases are fit 
using the Knudsen number formulation, and provide a 
reasonable fit to the theory points. We note that other 
parametcrizations described in give equally good fits. 
Then utilizing the information on the initial conditions 
in the model for R and assuming a temperature T = 200 
MeV (as was done in [i^), we extract the values for 77/s 
as a function of centrality (shown in the right panel of 
Figure HI). 

It it is immediately clear that the extracted values 
for 77/s are strongly centrality dependent, despite the 
fact that 77/s was input constant for all centrali- 

ties. While it has been postulated that this variation with 
centrality mimics viscous effects due to the finite lifetime 
of the system, this hypothesis then raises the question 
why such a large early freeze-out effect would maintain 
the expected curvature relation for the approach to ideal 
hydrodynamics as a function of centrality. Even more 
importantly, the numerical value for the 77/s extracted 
in central A-A collisions is incorrect by a factor of two. 
Taken together, these observations argue strongly against 
simply "calibrating" the procedure to remove the factor 
of two discrepancy. Any such renormalization should 
also be required to reproduce the input constancy of 
77/s, which in turn would require additional centrality- 
dependent adjustments to the parameters that appear in 
Equation [TJ Such an ad hoc procedure is clearly unsat- 
isfactory, especially given the availability of the very hy- 
drodynamic calculations that would be employed to reach 
a forced consistency. The inescapable conclusion is that 



the most reliable method currently available for extrac- 
tion of 77/s from experimental data is direct comparison 
to the output of hydrodynamic calculations, keeping in 
mind the issues identified in the preceding sections of this 
paper. It should be noted that there are additional con- 
cerns not addressed here, in particular that of transport 
in the hadronic phase. However, it is understood that 
the larger mean free paths in the hadronic phase lead to 
larger viscosity, so that the correct interpretation of the 
values used in the hydrodynamic calculations that ignore 
transport after hadronization would be upper bounds for 
77/s. 

VII. SUMMARY 

In this paper, we have explored some of the various 
inputs regarding hydrodynamic modeling of heavy ion 
collisions and the sensitivity of the final measured V2 
versus pt to these inputs. This work has been moti- 
vated by the striking agreement of V2 measurements over 
an enormous range in energy. Quantifying the precise 
level of agreement, and then in turn understanding the 
physics implications of the trends in the data requires a 
detailed understanding of these sensitivities. Even at the 
initial stage in such a systematic exploration, the argu- 
ment that the agreement is the result of fortuitous can- 
cellation of effects should be viewed with great scrutiny. 
As but one example, the agreement in the V2(j>t) data 
between RHIC and LHC is inconsistent with the predic- 
tions of models [30| which assume significantly different 
transport regimes at these two energies. It will be most 
interesting to extend hydrodynamic calculations to com- 
pare with new data at the lower colliding energies, al- 
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though issues of baryon contributions and modeling of 
the EOS will also need to be addressed. It is likely that 
these dynamic effects will greatly exceed the relatively 
modest variations in initial state eccentricity, which is no 
more than 10% when evaluated at the as a function of 
centrality rather than number of participants. 

We have also highlighted the need to take ratios which 
reveal that the greatest sensitivity to viscosity effects is at 
lower pt where the data is also the most accurate. These 
results highlight the need for more systematic studies 
and a re- evaluation of previously stated sensitivities to 
the early time dynamics and properties of the medium. 
Similarly, comparisons of the data to ever more sophisti- 
cated hydrodynamic and transport calculations, rather 
than parameterizations with underlying dynamical as- 
sumptions unsupported by detailed examination, are far 
more likely to lead to precision extraction of transport 
coefficients. An important development in these investi- 
gations is the public availability of modern and reliable 
hydrodynamic codes, which greatly leverage's the com- 



munity's ability to pursue these exciting topics. 



VIII. ACKNOWLEDGEMENTS 

We acknowledge useful discussions with Javier Al- 
bacete, Adrian Dumitru, Kevin Dusling, Matthew 
Luzum, Michael McCumber, Jean- Yves OUitrault, Paul 
Romatschke, and Raimond Snellings. We particularly 
appreciate the public availability of these algorithms and 
experimental data. JLN acknowledges funding support 
from the United States Department of Energy Division 
of Nuclear Physics Grant DE-FG02-00ER41152, the Dis- 
covery Center at the Niels Bohr Institute, and the Nordea 
Foundation. IGB is supported by the Danish National 
Science Research Council and the Danish National Re- 
search Foundation. WAZ is supported by U.S. Depart- 
ment of Energy grant DE-FG02-86ER40281. 



[22] 



K. Adcox et al. [ PHENIX Collaboration ], Nucl. Phys. 

A757, 184-283 (2005). [nucl-ex/0410003]. 

D. T. Son, A. O. Starinets, Ann. Rev. Nucl. Part. Sci. 

57, 95-118 (2007). [arXiv:0704.0240 [hep-th]]. 

Clifford V. Johnson and Peter Steinberg, Physics Today 

63, 5, 29-33 (2010). 

J.E. Thomas, Physics Today 63, 5, 34-37 (2010). 

B. Jacak and P. Steinberg, Physics Today 63, 5, 39-43 

(2010). 

B. Alver, G. Roland, Phys. Rev. C81, 054905 (2010). 
[arXiv: 1003.0194 [nucl-th]]. 

J. -Y. OUitrault, Phys. Rev. D46, 229-245 (1992). 

P. Sorensen, J. Phys. G G37, 094011 (2010). 

[arXiv: 1002.4878 [nucl-ex]]. 

K. Aamodt et al. [ The ALICE Collaboration ], 
[arXiv:1011.3914 [nucl-ex]]. 

N. Borghini, P. M. Dinh, J. -Y. OUitrault, Phys. Rev. 

C64, 054901 (2001). [nucl-th/0105040]. 

B. I. Abelev et al. [ STAR Collaboration ], Phys. Rev. 

C77, 054901 (2008). [arXiv:0801.3466 [nucl-ex]]. 

L. Kumar for the STAR Collaboration, arXiv:1101.4310 

[nucl-ex] . 

A. Adare et al. [ PHENIX Collaboration ], Phys. Rev. 
C77, 064907 (2008). [arXiv:0801.1665 [nucl-ex]]. 
R. Snelling, private communication. 

A. Adare et al. [PHENIX Collaboration], Phys. Rev. 
Lett. 105, 062301 (2010) [arXiv: 1003.5586 [nucl-ex]]. 
R. A. Lacey, A. Taranenko, N. N. Ajitanand et al., 
[arXiv:1011.6328 [nucl-ex]]. 

J. -Y. OUitrault, A. M. Poskanzer, S. A. Voloshin, Phys. 

Rev. C80, 014904 (2009). [arXiv:0904.2315 [nucl-ex]]. 

F. Cooper, G. Frye, Phys. Rev. DIO, 186 (1974). 

U. W. Heinz, [arXiv:0901.4355 [nucl-th]]. 

D. A. Teaney, [arXiv:0905.2433 [nucl-th]]. 

P. Romatschke, Int. J. Mod. Phys. E19, 1-53 (2010). 

[arXiv:0902.3663 [hep-ph]]. 

T. Hirano, N. van der Kolk, A. Bilandzic, Lect. Notes 
Phys. 785, 139-178 (2010). [arXiv:0808.2684 [nucl-th]]. 



[23] M. Luzum, P. Romatschke, Phys. Rev. C78, 034915 

(2008) . [arXiv:0804.4015 [nucl-th]]. 

[24] P. Romatschke, U. Romatschke, Phys. Rev. Lett. 99, 

172301 (2007). [arXiv:0706.1522 [nucl-th]]. 
[25] M. Luzum, P. Romatschke, Phys. Rev. Lett. 103, 262302 

(2009) . [arXiv:0901.4588 [nucl-th]]. 

[26] M. Laine and Y. Schroder, Phys. Rev. D 73 (2006) 

085009 [arXiv:hep-ph/0603048]. 
[27] C. Shen, U. Heinz, [arXiv:1101.3703 [nucl-th]]. 
[28] K. DusUng, G. D. Moore, D. Teaney, Phys. Rev. C81, 

034907 (2010). [arXiv:0909.0754 [nucl-th]]! 
[29] J. Jia, W. A. Horowitz, J. Liao, 

[arXiv:1101.0290 [nucl-th]]. 
[30] H. Niemi, G. S. Denicol, P. Huovinen et al., 

[arXiv:1101.2442 [nucl-th]]. 
[31] P. Huovinen, D. Molnar, Phys. Rev. C79, 014906 (2009). 

[arXiv:0808.0953 [nucl-th]]. 
[32] Y. Hidaka, R. D. Pisarski, Nucl. Phys. A820, 91C-94C 

(2009). [arXiv:0811.1202 [hep-ph]]. 
[33] H. B. Meyer, [arXiv:1012.0234 [hep-lat]]. 
[34] B. Alver, M. Baker, C. Loizides et al, [arXiv:0805.4411 

[nucl-ex]]. 

[35] B. Alver et al, Phys. Rev. C 77, 014906 (2008) 
[arXiv:071 1.3724 [nucl-ex]]. 

[36] S. A. Voloshin, A. M. Poskanzer, A. Tang and G. Wang, 
Phys. Lett. B 659, 537 (2008) [arXiv:0708.0800 [nucl-th]]. 

[37] T. Hirano, P. Huovinen, Y. Nara, [arXiv: 1012.3955 [nucl- 
th]]. 

[38] H. Song, S. ABass, U. WHeinz et al, [arXiv:1101.4638 
[nucl-th]]. 

[39] B. Schenke, S. Jeon, C. Gale, [arXiv: 1009.3244 [hep-ph]]. 
[40] H. - J. Drescher, A. Dumitru, A. Hayashigaki et al, Phys. 

Rev. C74, 044905 (2006). [nucl-th/0605012]. 
[41] J. L. Albacete, A. Dumitru, [arXiv:1011.5161 [hep-ph]]. 
[42] C. Alt et al [ NA49 CoUaboration ], Phys. Rev. C68, 

034903 (2003). [nucl-ex/0303001]. 
[43] H. -J. Drescher, A. Dumitru, C. Gombeaud et al, Phys. 

Rev. C76, 024905 (2007). [arXiv:0704.3553 [nucl-th]]. 



[44] J. L. Nagle, P. Steinberg, W. A. Zajc, Phys. Rev. C81, A830, 463C-466C (2009). [arXiv:0908.0403 [nucl- 

024901 (2010). [arXiv:0908.3684 [nucl-th]]. 
[45] H. Masui, J-Y. Ollitrault, R. Snellings et al, Nucl. Phys. 



0.8 



c 
« 

2" 0.7 
c 

(0 
Q. 

~ 0.6 

(0 
Q. 



cJ" 0.5 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



CO 



0.4 — 
0.3^ 
0.2 — 
0.1 




Monte Carlo Glauber 

^ — Pb-Pb2.76TeV (a^^ = 64 mb) 

Au-Au 200 GeV (a^^ = 42 mb) 
Au-Au 39 GeV (a^^ = 34 mb) 




J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L_ 



10 20 30 40 50 60 70 80 90 fOO 



Centrality Percentile (by N^^^^) 



O" 0.2 ^ — I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 





Impact Parameter (fm) 



